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Cosmology of Hofava-Lifshitz f(R) Gravity 
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We discuss the cosmology of recently proposed Hofava-Lifshitz f(R) gravity. In particular, we 
derive the modified Hubble equation that reduces to the standard HL gravity case in appropriate 
limit. We show how the bounce solutions in this theory are modified due to nonlinear effect of f(R) 
gravity. We also show that in principle, the Universe in this set-up can show cyclic behavior. 



I. INTRODUCTION 

Recently Hofava has proposed a theory of quantum 
gravity which is power-counting renormalizable in the 
ultraviolet (UV). This has been achieved through an 
anisotropic scaling between space and time, and there- 
fore it violates Lorentz invariance in the UV The in- 
frared limit of the theory reproduces General Relativity 
for a particular choice of a parameter, namely A = 1. The 
breaking of the Lorentz symmetry is performed by a pre- 
ferred foliation by three dimensional spacclike hypersur- 
faces, which in turn splits the coordinate into space and 
time. It allows one to write the Einstein Hilbert action of 
the General Relativity with higher spatial derivatives of 
the metric. This improves the UV behaviour of the gravi- 
ton propagator and makes the theory power-counting 
renormalizable. Moreover, the action has only second 
order time derivatives preventing the presence of ghosts 
in the theory. However, the explicit breaking of the gen- 
eral covariance by such foliation of the space-time intro- 
duces a new scalar degree of freedom which shows some 
unavoidable pathological behaviours Q. Subsequently, 
there were many more proposals to improve the theory 
and overcome these pathological behaviours Q . It is also 
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worth mentioning here, that at the one loop quantum 
level-the violation of Lorentz invariance has dramatic ef- 
fect 0, Q and in order to have a viable Hofave-Lifshitz 
(HL) gravity much works need to be done towards un- 
derstanding the renormalization group flows. 

Apart from the abovementioned problems in the HL 
theory, it still is interesting to look at different aspects of 
this theory from the perspective of cosmology. Many in- 
teresting cosmological implications have been found for 
the HL gravity |6|. It was found out that Friedmann 
equations in HL gravity contains some terms which be- 
haves as dark radiation with negative coefficient. There- 
fore it is possible to find bounce solution which has a 
non-singular behaviour of scale factor as the origin of the 
Universe • Different UV scaling of HL gravity can also 
lead to scale-invarint cosmological perturbations without 
any need of inflation Q . Given the fact that HL gravity is 
a non-rclativistic modified gravity, the idea of cyclic cos- 
mology was first realized in this non-relativistic gravity 
by Cai and Saridakis Q . In this work, only the gradient 
terms of the gravitational action were modified and the 
model was free from the strong coupling problem which 
was extensively discussed in Horava-Lifshitz models [l[. 

On the other hand, many recent observations indicate 
that the Universe has probably gone through two phases 
of accelerated expansion: one in the early epoch to solve 
many cosmological puzzles with the creation of observed 
density perturbations (lfjj |. another to accommodate the 
peculiar observations that the current Universe is accel- 
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erating ll| . Instead of adding some extra matter compo- 
nent in the Universe to invoke the accelerated phases, an 
alternative is the modification of Einstein-Hilbert term 
by replacing the Ricci scalar R in the action with some 
general functions f(R) of the Ricci scalar. For an ex- 



tensive review see [12j. Considering that HL gravity is 
a promising candidate for an UV completion of General 
Relati vity , / (it!) versions of HL gravity were also pro- 



posed [13l - tl5( | . 

In this paper we make an analysis of the cosmology of 
f(R) HL gravity that has been proposed in (lij . In par- 
ticular, we will derive the modified Hubble equation that 
reduces to the standard HL gravity case in appropriate 
limit. Due to nonlinear effect of f(R) gravity we will see 
how the bounce solutions are modified. We will explicitly 
show that in early time the Universe can go through a 
bounce phase, whereas in late time it admits turn-around 
solution. Therefore, in principle, the Universe in this set- 
up can show cyclic behavior. After describing the basics 
of HL gravity and its cosmology in Sec. 2 and Sec. 3 
respectively, we discuss a particular form of f(R) HL 
gravity in Sec 4. In Sec 5 we discuss the cosmology of 
this particular scenario. 



II. A BRIEF REVIEW OF HORAVA-LIFSHITZ 
GRAVITY 

In this section we briefly review the Hofava-Lifshitz 
(HL) gravity as was originally proposed in [l[. The dy- 
namical variables in HL gravity are similar to the ADM 
formalism of Einstein gravity, namely the lapse function 
N, shift vector JVj and the spatial metric g^. We con- 
sider only the projcctablc version of HL gravity where 
the lapse function needs to be function of time t only. In 
terms of these variables the metric is 



ds 1 



N 2 dt 2 + g tJ {dx % + N l dt){dx J + N 3 dt), (1) 



where Ni = gijN^ . The gravitational part of the action 
consists of a kinetic and a potential part 



S g = J dtd 3 xy^N (C K + Cv) ■ 



The kinetic terms are given by 



C K = -^(KijKV -\K 2 ), 



(2) 



(3) 



with covariant derivatives defined with respect to the 
spatial metric gij. The potential term in the "detailed- 
balance" 1 form is given by 
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(5) 



where the super-metric Gij-ki depends on an arbitrary di- 
mensionless coupling constant A (susceptible to quantum 
corrections) 

Qij-M = ^{9ik9jl + 9il9jk) + Y~^9ij9kh (6) 
The E tensors are given by 

E ij = 2^ y y) (7) 

w z 2 



where 



(8) 



is the Cotton tensor, a conserved traceless tensor that 
vanishes for conformally flat metrics, k, w, [i are con- 
stants with mass dimension —1, 0, 1 respectively. A^ is 
related to the cosmological constant in the IR limit of 
the theory. 

Demanding that in the infrared limit, the HL action 
reduces to the Einstein-Hilbert action 



S 



EH 



ike J <''- r \ 7 l Xih -J K ' ! -K 2 + R- 2 A) (9) 



allows us to write the speed of light, Newton's constant 
and the cosmological constant in terms of parameters of 
the model 



K [l / A 



G=4~S A = | A - ( 10 ) 
6Zirc 2 



4 V 1 -3A' 



In addition we have to identify A = 1. From Eq.(|10p we 
can easily see that for A > i, the cosmological constant 
A w must be negative. 



III. COSMOLOGY OF HL GRAVITY 



Now we discuss the implications of HL gravity in the 
cosmological context 0,00- In particular, we assume 



where the extrinsic curvature is given by 
1 



K 



13 ~ 2N^ ij ~ ViNj ~ VjNi) ' 



, . 1 Without the detailed balance condition, the most general renor- 
l J malizable theory contains cubic order terms in curvature F^. ll6H . 
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the homogeneous and isotropic Universe described by the 
following FRW metric 



gij dx l dx^ = a(ty 



dr 2 



1 — kr 2 



sin 9dcj) 



(11) 



with N = 1, N l = 0, where k = +1,0,-1 corresponds 
to a closed, flat or a open Universe respectively. Due 
to many symmetries of the FRW space-time the calcula- 
tions simplify by a large amount, but still leaving some 
novel signatures of the HL gravity. Assuming matter con- 
tribution behaves as an ideal fluid we find the following 
Hubble equation by varying the action with respect to 
the lapse function N 



k 





2A w a 4 



(12) 



Variations with respect to gives us the following equa- 
tion 



H + -H 2 = —L 

2 3A - 1 V 2 



^ ~ ~ 2A^ ' ' (13) 



where p is the pressure of the fluid. 

For the vacuum case with p = p = 0, and we have 



H 2 = 



A, 



3A- 1 



1 - 



A,„a 2 



(14) 



As we have noted earlier, for A > i, the bare cosmo- 
logical constant A w that appears in the HL theory must 
be negative. Therefore the only solution (static) exists 
for k = — 1 with a 2 = — 1/ A w . However one can can also 
make an analytic continuation to extend the cosmological 
constant into positive regime in the frame of HL gravity. 
This issue has been already studied explicitly in Q. In 
our present investigation, we will only consider the case 
where the bare cosmologial constant is negative. 

We find modifications to the standard Hubble equation 
due to higher order derivative terms in the action. Inter- 
estingly enough, these contributions from higher deriva- 
tive terms vanish for k = 0. In addition, they are im- 
portant only at small a, giving us back to the standard 
cosmological solutions in the IR. We also see that the 
Hubble equation is modified by an overall factor depend- 
ing on A. 

The most important modification is the term propor- 
tional to 1/a 4 , that contributes a negative energy density 
for any nonzero spatial curvature. This term is similar 



to the "dark radiation" in the braneworld scenario, with 
an important difference of being nonzero only for nonzero 
spatial curvature [171 ] . If in a collapsing phase (decreasing 
a) of the Universe at a very early time, the matter energy 
density grows slower than a -4 , the dynamics leads to a 
bounce with H = 0, and H > [ljj ]. The condition for 
bounce in HL gravity turns out to be 2 . 



(I--) 



-V >0. 



(15) 



During the contracting phase of the Universe, matter can 
be described by an oscillating scalar field with a potential 



1 



V(<j>) = -m 2 (j) 2 



(16) 



During the contraction phase, the amplitude of oscilla- 
tion grows until it reaches a critical value where energy 
density carried by the field becomes constant. Subse- 
quently, due to the rapid contraction of the Universe, 
negative contributions coming from 1/a 4 term catches 
the energy density of <f> field. In the end Universe passes 
through a bounce when Hubble constant H = with 
H > 0. The application of this mechanism to achieve a 
nonsingular bounce was firstly appeared in the model of 
Lee Wick theory 0. 

In the following section we will find out the cosmolog- 
ical implications of modified HL gravity. 



IV. f(R) HORAVA-LIFSHITZ GRAVITY 

As an explanation of early and late time accelera- 
tion of the Universe, modifications of Einstein gravity 
without an extra matter source is a viable option (for 
a general review see [l^l)- In this approach, the stan- 
dard Einstein-Hilbert action is added/replaced by the 
higher order curvature invariants such as R 2 1 R lLU W ,v 
etc. Indeed, there are examples of f(R) gravity theo- 
ries where the early time inflation can be unified with 
late time cosmic acceleration [iof . On the other hand, 
Hofava-Lifshitz Gravity is a promising way forward in 
constructing the theory of quantum gravity 0. The the- 
ory is power-counting renormalizable in the expense of 
breaking Lorentz symmetry explicitly. The ultraviolet 
behavior of the theory is improved by the higher spatial 
derivative terms of the metric. Subsequently, f(R) HL 
gravity has been constructed recently following the way 



2 Condition is different than what has been found in [IS 



4 



the original Einstein gravity was modified to f(R) theo- 
ries. In this section we review these approaches briefly. 

A straightforward modification of HL gravity has been 
recently proposed by 0, [lfi] 



S = / dtd 3 x^gNf(R), 



(17) 



with 



R 



K^Kij - XK 2 + 2fiV ' v n v ~ n v V v rt 



-E^g im E 



kl 



(18) 



where is a unit vector perpendicular to i = con- 
stant space-like hypersurface T, t . Similar to the stan- 
dard Einstein-Hilbert action, the term proportional to p, 
is usually dropped in HL gravity as it is a total derivative 
term in the action. But in the non-linear generalization 
to f(R) gravity, this term can not be dropped (l5j. Mod- 
ified version of this HL gravity has been subsequently dis- 
cussed for unifying inflation with dark energy 



21 



15 



and 



22|. 



Hamiltonian analysis has been carried out in 

Following Hofava's proposal, a different version of 
modified f(R) HL gravity has been introduced jl3| early 
on. In this case, at first a partition function of D- 
dimensional gravity theory was constructed. Demand- 
ing the existence of a (D + 1) dimensional quantum 
gravity theory such that the norm of its ground state 
wave function coincides with the partition function of D- 
dimcnsional theory, it was then shown that infinite num- 
ber of Hamiltonian obeys the detailed balance conditions. 
For a special form of the Hamiltonian which corresponds 
to f(R) = V~R, a Lagrangian formulation of the theory 
was also made assuming projectability condition. Action 
for this modified HL gravity is given by [13 1 



S = 



dtd 3 




M 6 ziKyKV - XK 2 ) - M 6 (19) 



Here M is a mass parameter. In this form of the ac- 
tion, the dimensions of different quantities are as follows: 



[G 



ijkl\ 



[Mass] ,^] 



[Mass] 3 , [k 2 ] 



[Mc 



and 



[Kij] = [Mass] 1 . In the original paper by Kluson 131 ] . 
the action was written in a unit where M = 1. In our 
subsequent calculations, we also assume M = 1. 

A similar form of the action was also found without 
projectability condition. In a subsequent paper a gener- 
alization for any form of f(R) was presented [14j. The 
new form of f(R) HL gravity are invariant under folia- 
tion preserving diffeomorphism similar to the original HL 



gravity 0. We note that the ordinary f(R) gravity the- 
ories have full diffeomorphism invariance. It is trivial to 
check that linearized version of the above action in Eq. 
(|19p gives us back the original HL action. Goal of this 
article is to present the cosmological consequences of the 
modified f(R) HL action given by Eq. (fT9l) . 



V. COSMOLOGY OF f(R) HORAVA-LIFSHITZ 
GRAVITY 

In this section we will discuss the cosmological solu- 
tions of modified f(R) HL gravity defined by Eq. (JT5J. 
In particular, we will look for the 'bounce' and 'turn- 
around' solution in this scenario and point out the main 
differences with the standard HL gravity solutions. 

Assuming the background Universe described by a 
FRW metric and the matter contribution is equivalent 
to an ideal fluid described by energy density p and pres- 
sure p, the Hubble constant can be extracted as 



1 



6(3A- 1) 



1 - 



K 2 p 



El 

2 



6A„ 



12- 



6 



A™ a 

(20; 



This is the modified Hubble equation for f(R) HL grav- 
ity and should be compared with Eq. (|12|) . We note that 
the vacuum solution with energy density p = is exactly 
similar to the standard HL gravity theory. Therefore, 
similar to the standard HL gravity, vacuum solution ex- 
ists only for open Universe (k = —1) with scale factor 
a = y/—l/A w , where A w is negative for A > 1/3. Thus, 
this particular form of f(R) gravity does not change the 
pure gravity part in FRW space-time. In contrast to the 
other version of f(R) HL gravity [15j, the vacuum solu- 
tion does not admit de-Sitter solution. As we see from 
Eq. (|19[) , similar to the standard HL gravity, many novel 
terms are proportional to the spatial curvature. Another 
important feature of the modified Hubble equation is that 
even for a flat case (k = 0), there is a modification in the 
Einstein equation. This is because of the presence of p 2 
term. Moreover this term comes with a negative sign. 

Expanding the above expression for the Hubble con- 
stant in terms of p/k 2 we find 



H 2 



2 (K, 
3A - 1 V 2 

3 A w p 2 
2 K 4 

3p 2 k 2 
+ 2a 4 K 4 A,, 



k 2 



A„ 



' 12' 
2kp 



2A,, 



3kp 2 



k 2 p 



a 4 K 2 A„ 



0( P 3 ) 



(21) 
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Firstly, in the leading order we recover the Hubble equa- 
tion for standard HL gravity. Higher order terms are due 
to nonlinear f(R) HL gravity, and many of them con- 
tributes to the "dark radiations", i.e redshifts as a -4 . 

For a flat universe with k = 0, from the positivity of 
H 2 we can easily find out that the modified theory must 
satisfy the following condition 



k 2 - y^K 2 - 12|AJ < p < k 2 + yjn 2 - 12|AJ (22) 

with |Au,| < k 4 /12. Thus the energy density p coming 
from the matter sector is both bounded above and below. 
This condition is starkly different than the standard HL 
gravity where we find a lower limit on the energy density 
6 1 A™ | < K 2 p for a given value of A w . 

Note that the Hubble parameter blows up at the singu- 
larity p — k . Although H is also singular at this point, 
the scale factor, energy density and pressure remains fi- 
nite. This is a new _kind of singularity and can not be 
classified following 



231. But we should stress that this 



singularity occurs at the Planck's scale where the classi- 
cal FRW equation for the cosmological evolution should 
not be valid and one should consider the quantum grav- 
ity effects. So for our purpose, working in the classical 
regime p < k 2 , avoids hitting this singularity. We com- 
ment on the avoidance of this singularity while having 
bounce solution in the next section. 

In the following we will discuss both the early and late 
time cosmology that is realized with this particular form 
of f(R) HL gravity. In particular, we will see, whereas in 
early time the Universe can go through a bounce phase, 
in late time it admits turn-around solution. In summary, 
the Universe can go through a cyclic behavior subject to 
some conditions satisfied by the matter content of the 
Universe. 



A. Bounce Solution 

One of the important aspect of HL gravity is that it 
allows the Universe to have bounce solutions without in- 
troducing any exotic matter components. As we have 
discussed earlier, this happens naturally due to "dark ra- 
diation" term in the Hubble equations. In the modified 
f(R) HL theory bounce structure changes mainly due to 
the presence of an extra term in the Hubble equation 
that depends on the square of energy density and has a 
negative coefficient. Now, in addition to the standard en- 
ergy density term, this term also contributes in making 
the Hubble constant zero as Universe passes through a 
contraction phase. 



Varying the action with respect to <?y gives us the ex- 
pression for the time variation of Hubble constant as 



H 



1 



12(3A- 1) 



3( P + P)(P- 

-2 'Kp + p) 



H 



24fc 
"a 5 " 



2Ak 2 
\A w \a 

-l 



(23) 



It is important to remember that A w is a negative quan- 
tity. A bounce is defined by the condition of H = 0, with 
H > 0. Around the bounce the term proportional to A^ 
is irrelevant, whereas the term proportional to a -2 grows 
less faster than a -4 term, thus we neglect those. As we 
have noted earlier, the bounce solution is possible only 
when the Universe is not spatially fiat (k ^ 0), and the 
solution remains same for both open or close Universe. 
Setting H = therefore allows us to find the constraint 



pi 

2 



6k 2 



o. 



(24) 



where p* and a* are the energy density and the scale fac- 
tor respectively when bounce happens. Now, demanding 
the positivity of H gives us the bounce condition as 



p: 



K 2 p* 



+ 3p*(p* — k ) > 0. 



(25) 



If we describe the matter with equation of state param- 
eter w with pressure p = wp, the above condition can be 
re-casted conveniently as the following 

p4l + 3w) + K 2 (l-'iw) >0. (26) 

Note that (f2"4"| has two solutions pf = n 2 ± 
pr. If we assume the negative sign for the 



critical density then the bounce occurs at the value of the 
energy density which is always smaller than the value at 
the singularity p = k 2 mentioned in the previous section. 

As we can see, a large range of values of w can easily 
satisfy the bounce condition. The novel feature of having 
terms in the Hubble equation proportional to a saves 
us from introducing exotic matter for realizing bounce 
in HL gravity scenario. The most important difference 
between the condition of Eq. (|26l) and the condition for 
standard HL gravity bounce is the appearance of higher 
order terms as p 2 . If we model the bounce by a scalar 
field with nearly flat potential, we have w ~ — 1, and in 
that case the condition can be written in the simple form 
p* < 2k 2 . For this case also, we can always avoid the 
singularity at p = n 2 mentioned earlier but can have the 
bounce, if we ensure that the evolution begins after the 
singularity i.e. p < k 2 . For w = and w = 1/3, one 
can similarly check that the singularity at p = k 2 can be 
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avoided. As we have seen earlier, for the standard HL 
gravity case, a bounce is always guaranteed if p ~ —p, 
but for the modified HL case, there is an upper bound 
on the total energy density of the bouncing matter field. 

After finding the conditions for the existence of bounce 
solutions, we now numerically solve the Hubble equa- 
tion to show the dynamical behavior of a bounce. The 
plots for the scale factor and the Hubble constant are 
shown in Fig. ([I]) and ([2]) respectively when energy den- 
sity is assumed to be a cosmological constant. Without 
the existence of any exotic matter content, the Universe 
shrinks to a minimal value of the scale factor and then ex- 
pands again. At the bouncing point, the Universe evolves 
through zero Hubble constant. 



0.52 - 



0.48 - 



l).4(i 



FIG. 1: The scale factor a(t) in the bouncing solution (for 
w = — 1) shrinks to a non zero value at the bouncing point 
and expands again. 




B. Turnaround Solution 

As we just discussed in the previous subsection, HL 
gravity naturally admits bounce solution if the Universe 
starts to contract at some point. In a similar way, the 
structure of the HL gravity is such that in its expand- 
ing cosmological phase, the Universe does not expand 
forever. For large values of the scale factor the Universe 
make a turnaround with Hubble constant being zero with 
H < 0. 

Turnaround solution can be easily seen in the standard 
HL gravity solution. For large value of the scale factor 
we can easily neglect all terms proportional to the spatial 
curvature, therefore 



H 



3A - 1 



I A, 



12 P 



(27) 



Now it is easy to see that turnaround always happens for 
non-zero value of cosmological constant with equation of 
state parameter w > — 1 at = ■ 

In a similar fashion, in the limit of large scale factor, 
the Hubble constant for the particular f(R) HL gravity 
that we are considering can be written as 



E- 



1 



6(3A- 1) 



K 2 p ■ 



6A„ 



f - 



(28) 



Note that this is valid only for the flat universe with k = 
0. For non flat universe, we can safely ignore p 2 terms 
in late times for getting turn around solutions which is 
similar to the standard HL scenario. In this case, the 
condition of H = at late time gives us following two 
solutions for critical turnaround density 



H(t) 

4 - 




FIG. 2: The Hubble parameter H(t) in a bouncing scenario 
(for w = — 1) showing that it crosses zero when the universe 
goes from a contracting phase to the expanding one. 



i± \ i 



12|A t 



(29) 



The condition for turnaround with H < can be written 
in general as the following 



(1 +w)(p* -k 2 ) <0 



(30) 



and can be satisfied with either w < — 1 and > k 2 , 
or w > —1 and /?* < k 2 . It is clear that the first condi- 
tion corresponds to the solution pj, whereas the second 
condition corresponds to p~ . Firstly, for A w = the 
turnaround density becomes p+ = 2k 2 , whereas other 
solution becomes trivial. In this case, the condition for 
bounce is w < — 1, i.e the Universe must be dominated 
by some phantom like field. We note that there is no 
turnaround solution for standard HL gravity case with 
A u , = 0. In our case, this feature is due to the presence 
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of negative energy density square term in the Hubble 
equation. For A w 7^ 0, the relevant solution is p~ with 
the condition for bounce being w > — 1 and p* < k 2 , and 
it can be achieved by any matter field satisfying strong 
energy condition. 

In summary, in this version of modified f (R) HL gravity 
the Universe can go through alternate cycle of contrac- 
tion and expansion via bounce and turnaround solution. 
This is possible even if the cosmological constant A w van- 
ishes in our model. 

VI. CONCLUSIONS 

In conclusion, we have studied the cosmology for the 
f(R) gravity model proposed recently by Kluson fljj ). 
Our main motivation is to sec whether one can get bounc- 
ing as well turnaround solution in this model. There 
are several interesting features in this model. Firstly the 
cosmology in the pure gravity case is exactly similar to 
what one gets in a standard HL gravity. But in the pres- 
ence of matter, things change substantially. There is a 
new correction term proportional to p 2 which also comes 
with a negative sign. Due to the presence of such term, 
the cosmology is different from the standard one even in 
the flat case. This is unlike the pure HL gravity, where 
modification to the standard cosmology only happens 



for a nonzero spatial curvature. Both the bounce and 
turnaround conditions in our model are different from 
pure HL gravity. Specifically, we need no exotic fluid to 
get the bounce; one can achieve this with normal mat- 
ter like radiation. Similarly turnaround can be achieved 
even in the absence of a cosmological constant although 
one may need a phantom fluid for such behaviour. 

Overall the f(R) HL gravity has many interesting cos- 
mological signatures as far as the background evolution 
of the universe is concerned. It will be interesting to see 
how the inhomogeneities grow in this model and this will 
be addressed in future works. 
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